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Context. Accurate estimates of the errors on the cosmological parameters inferred from cosmic shear surveys require accurate esti- 
mates of the covariance of the cosmic shear correlation functions. 

Aims. We seek approximations to the cosmic shear covariance that are as easy to use as the common approximations based on normal 
(Gaussian) statistics, but yield more accurate covariance matrices and parameter errors. 

Methods. We derive expressions for the cosmic shear covariance under the assumption that the underlying convergence field follows 
log-normal statistics. We also derive a simplified version of this log-normal approximation by only retaining the most important terms 
beyond normal statistics. We use numerical simulations of weak lensing to study how well the normal, log-normal, and simplified 
log-normal approximations as well as empirical corrections to the normal approximation proposed in the literature reproduce shear 
covariances for cosmic shear surveys. We also investigate the resulting confidence regions for cosmological parameters inferred from 
such surveys. 

Results. We find that the normal approximation substantially underestimates the cosmic shear covariances and the inferred parameter 
confidence regions, in particular for surveys with small fields of view and large galaxy densities, but also for very wide surveys. In 
contrast, the log-normal approximation yields more realistic covariances and confidence regions, but also requires evaluating slightly 
more complicated expressions. However, the simplified log-normal approximation, although as simple as the normal approximation, 
yields confidence regions that are almost as accurate as those obtained from the log-normal approximation. The empirical corrections 
to the normal approximation do not yield more accurate covariances and confidence regions than the (simplified) log-normal approx- 
imation. Moreover, they fail to produce positive-semidefinite data covariance matrices in certain cases, rendering them unusable for 
parameter estimation. 

Conclusions. The log-normal or simplified log-normal approximation should be used in favour of the normal approximation for pa- 
rameter estimation and parameter error forecasts. More generally, any approximation to the cosmic shear covariance should ensure a 
positive-(semi)definite data covariance matrix. 
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1. Introduction 

In recent years, observations of weak lensing by large-scale 
structure, called cosmic shear, have become an important tool 
for studying the Universe. Remarkable constraints on the mat- 
ter content of the Universe, its expansion history, and the am- 
plitude of the cosmic density fluctuations have been obtained, 
e.g., from measurements of cosmic shear in the Canada-France- 
Hawaii Telescope (CFHT) Legacy Surve y ([Semboloni et al. 



> |2006||Hoekstra et al.|2006||Fu et al |2008| , the lOOdeg^ weak- 



lensing survey ( Benjamin et al. 2007|l, a ndthe Cosmic Evolution 
Survey (COSMOS HMassey et al.|2007l|Schrabback et al.|20"T0l i. 

Future cosmic shear surveys are expected to provide essen- 
tial information on the properties of the dark matter and the 
dark energy, and on possible deviations of gravity from General 
Relativity (Huterer 2010). 

cosmic shear surveys aim at inferring the statistical proper- 
ties of the convergence and shear field from the observed im- 
age shapes of distant galaxies. There are various statistical mea- 
sures that probe the two-point statistics of the convergence and 
shear, e.g. the shear correlation functions (Blandford et al. 1991 ; 
[Kaiser 1992), the shear power spectra ( [Kaiser, 1992), the shear 
dispersion ( ^Mould et al. 1994 1, the aperture-mass dispersion 



( Schneider et al. 1998| l, the ring st atistics ([Schnei de r & Kilbinger 
'2007[), the statistics proposed by [Fu & Kilb inger (2010[), or the 
COSEBIs ( [Schneider et al.|2010 1. In practice, however, only the 
shear correlation functions ^+ and ^_ are estimated directly from 
the data (since they can be estimated most easily), and any other 
two-point statistics are computed from these estimates. Thus, ac- 
curate knowledge of the noise properties of the shear correlation 
estimators is needed for reliable estimates of the errors on the 
cosmological parameters inferred from cosmic shear surveys. 

The number of independent measurements in a cosmic shear 
survey is usually too small to permit an good estimate of the 
full covariance of the measured statistics directly from the data. 
Therefore, one has to resort to covariance estimates obtained by 
other means. Supposedly the most accurate, but also computa- 
tionally by far the most expensive approach is the direct esti- 
mation of a cosmic shear covariances from A^-body simulations 
of cosmic structure formation (e.g. Sato et al. 2011b, or this 
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work). This method requires many independent realisations of 
the survey field. So far, it has been successfully applied to lens- 
ing surveys covering a small fraction of the sky, e.g. the Chandra 
Deep Fie ld South (Hartlap et al.|[2009]) aiid the COSMOS field 
(Schrabback et al.i2010 ; Semboloni et al.|2011| l. However, cre- 
ating sufficiently many realistic mock observations to calculate 
covariances for future all-sky lensing surveys in various cos- 
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mologies and parameter ranges will remain a very challenging 
computational task in the foreseeable future. 

Since reliable estimates for the covariances of cosmic shear 
two-point statistics from real or mock data are often out of 
reach, it appears desirable to have simple approximations at 
hand, which are easy to compute and yield sufficiently accu- 
rate covariance matrices and parameter errors. General expres- 
sions for the covariance of lensing two-point statistics, including 
the shear correlation functions, have been derived by [Schneider] 
et al. ( 2002a i and'Joachimi et al. ( 2008 ) under the assumption of 
normal (Gaussian) statistics for the shear and convergence field. 
This normal approximation is relatively easy to apply, and re- 
quires no additional information besides the convergence corre- 
lation itself and the basic survey parameters. Simulations show, 
however, that the normal approximation tends to underestimate 
the covariance of two-point lensing statistics (e.g. Scoccimarro, 



etal. 


1999 


Cooray & Hu|2001 Van Waerbeke et al.|2001 Sato 


et al. 


2009 





There have been various attempts to improve the normal ap- 
proximation by amending the resulting expressions for the cos- 
mic shear covariance by scale-dependent correction factors that 
are calibrated with mock data from structure formation simula- 
tions (Van Waerbeke et al.|pOOT] [Semboloni eFaLll^OOTl [Sat^ 
|et al. 2011b) . A drawback of this approach is that, while any 
corrections to the normal approximation depend in a non-trivial 
way on the assumed cosmology, the correction factors can only 
be tuned to the particular cosmology of the simulations at hand. 

Other approaches towards more accurate lensing statistics 
are based on halo raodels jCooray & Hu |2001 I Takada & Jain 



[2009) [Pielorz et al.|[20T0l|Kainulainen & Marra|201 l| l. The re- 

sulting expressions for the covariance of lensing two-point statis- 
tics require considerable more effort to compute than the corre- 
sponding expressions in the normal approximation, even if sim- 
plified by fitting formulas ( |Pielorz et al.|2010) . A thorough as- 
sessment of the accuracy of the halo model approach to cosmic 
shear covariances (e.g. by comparing its results with those from 
large high-resolution simulations of weak lensing) is still out- 
standing. 

Normal random fields are arguably the most simple statisti- 
cal models for the three dimensional matter distribution and the 
resulting weak-lensing convergence and shear fields. However, 
as mentioned above, they fail to accurately describe those fields' 
higher-order correlations such as the covariance of two-point 
correlations. There are few other classes of random fields that 
allow one to calculate their higher-order correlations analyti- 
cally. These include log-normal random fields, on which we fo- 
cus here. 

Already p^ubble| p 934| l noted that observed galaxy densities 
roughly follow a log-normal distribution. Coles & Jones ( 1991| l 



suggested to describe the three-dimensional matter density dis 
tribution in our Universe by a log-normal random field. Later, 
numerical simulations of cosmic structure formation showed 
that a log-normal field provides indeed a good description of 
the one- and two-point statistics of the evolved matter density 
field, even though a point-wise exponential mapping of the initial 
Gaussian density field does not describe the evolved density field 
well (Kofman et al. 1994; K ayo et al.|2001[ ). The one-point prob- 
ability density function (pdf) of the projected density and lens- 
ing convergence are also fitted well by a log-normal pdf (|Taruya 



|et al.|2002l l. Even better fits to the convergence pdf, particularly 
in the tails of the distribution, can be obtained by generalisations 



In this work, we derive expressions for the cosmic shear co- 
variance under the assumption that the underlying convergence 
field follows log-normal statistics. Furthermore, we derive a sim- 
plified version of this log-normal approximation by only retain- 
ing the most important terms beyond normal statistics. In partic- 
ular the simplified log-normal approximation meets our aim of 
an expression for the cosmic shear covariance that is as simple 
as the normal approximation. 

We use numerical simulations to show that the log-normal 
and the simplified log-normal approximations reproduce the 
shear covariances for cosmic shear surveys much more accu- 
rately than the normal approximation. Furthermore, we demon- 
strate that, in contrast to the normal approximation, the uncer- 
tainties in the cosmological parameters inferred from the log- 
normal and the simplified log-normal approximation are in good 
agreement with the errors resulting from the covariances we es- 
timated from the numerical simulations. We also show that the 
covariances and parameter errors in the log-normal approxima- 
tion have accuracies comparable to or better than those derived 
from empirical corrections to the normal approximation pro- 
posed by Semboloni et al. ( 2007 1 and Sato et al. (201 lb ). In ad- 



dition, we discuss the problem that approximations to the cosmic 
shear covariance, in particular those based on empirical fits, may 
produce covariance matrices that are not positive-semidefinite, 
which makes their use in parameter estimation at least question- 
able. 

The paper is organised as follows. The theory of estimators 
for cosmic shear two-point statistics and their noise properties 
are briefly discussed in Sect. [2[ In Sect. [3] we describe our nu- 
merical simulations for creating mock lensing surveys and for 
studying the statistical properties of the convergence and shear. 
The results of the simulations and the comparison between the 
different approximations to the cosmic shear covariance are pre- 
sented in Sect. [4] The main part of the paper concludes with 
a summary and discussion in Sect. [5] A more detailed discus- 
sion about normal and log-normal random fields, and the covari- 
ance of estimators for the cosmic shear 2-point correlation can 
be found in the Appendix. 



2. Theory 

Here, we briefly discuss weak gravitational lensing, estimators 
for the cosmic shear correlation functions and their noise prop- 
erties, and the estimation of cosmological parameters from cos- 
mic shear data. The discussion of the statistical properties of the 
cosmic shear estimators is based on the real-space approach by 
Sch neider et al.| (|2002a), but includes log-normal convergence 
fields in addition to normal convergence fields. 



2.1. Gravitational lensing 

Gravitational lensing, the deflection of photons from distant 
sources by the gravity of intervening matter structures (e.g. 
Schneider et al. 2006) , causes shifts in the observed image posi- 
tions relative to the sources' 'true' sky positions. The shifts may 
be described by a deflection field 



(1) 



of a log-normal pdf ( 


Das & Ostriker 2006 Takahashi et al. 201 1 


IJoachimi et al.|20I I 





which relates the observed image position - (01,62) of a 
point-like source at redshift z to its true angular position /3 = 
(J3i{9,z),/32(0,z)). The distortions induced by diff'erential deflec- 
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tion can be quantified by the distortion matrix 



K + 7i 72 + w 
72 - w K-Jl 



(2) 



which can be decomposed into a convergence k(0, z), a complex 
shear y{0,z) - yi{0,z) + vy2{0,z), and an asymmetry io{0,z)- 

We define the effective convergence and shear for a popula- 
tion of sources with normalised redshift distribution p^iz) by 



k(0)^ jdzpz(z)K(e,z) and 
y(0)^ jdzp,Sz)j(0,z). 



(3) 
(4) 



The tangential component y^{0,-9) and cross component 
yx(0,tf') of the shear y(0) at position with respec t to the di- 
rection are defined by (e.g. Schneider et al.|2002b i 



yriO,^) = -1R (r(6')e-^'^<*^) and 



(5) 
(6) 



where (p{&) denotes the polar angle of the vector if. 

For a statistically homogeneous and isotropic universe, the 
convergence correlation function can be defined by the ex- 
pectation (e.g. |Blandford et al.|1991[|Kaiser|1992) i 



um^(K(0)K(0+^)). 



(7a) 



Here, {/) denotes the ensemble average of a function / over 
all realisations of the observable universe for a fixed set of cos- 
mological parameters. The two shear correlation functions ^+(j?) 
and ^-()?) are defined by 

^±{m) = {yT(0, »)7t(0 + ^, ^)} ± {yx(0, ^)yx(0 + ^, ^)} ■ (7b) 

If one neglects higher-order lensing effects (which are small, 
cf. [Hilbert eF al. 2009 ; Krause & Hirata 2010) and contamination 
by observational systematics, the deflection field Q is a gradi- 
ent field, the asymmetry co vanishes, and the convergence k and 
shear y are related by ( Kaiser|19^ i 



As a consequence ( Crittenden et al.|2002| l, 

= U-») and 



where 



(8) 

(9) 
(10) 



G-()?,i?') 



4 



127?'2\ , 1 



(11) 



Here, H denotes the Heaviside step function, and do denotes the 
Dirac delta function. 



2.2. Estimators for the shear correlations and their 
covariance 

A suitable estimator for the shear y{0^'^) at the sky position 0*'' of 
a galaxy / at redshift z*'* from a sample with redshift distribution 
p~{z) is the observed image ellipticity y^'K To lowest order, the 
observed ellipticity is a sum. 



i(0 Ji}\ 



(12) 



of the 'true' shear y(^*'\z*'') and an intrinsic ellipticity e*'\ 
which vanishes on average. From a survey with A^g galaxies pro- 
viding shear estimates -y*'', . . . at positions 0^^\...,0^^\ 
one can estimate the shear correlations ^±(i?) at separations § 



J- 



Tx 7 x 



Here, the bin window function 



A(i?, e) = 



1 for \6--&\< A(&)/2 and 
otherwise. 



(13) 



(14) 



where A{§) denotes the bin width, and ''^ and 7x"'' denote the 
tangential and cross component, respectively, of the shear esti- 
mated from the shape of galaxy / with respect to the line joining 
galaxy / and j. 

If (i) the bin width is small compared to the scales on which 
correlations change, (ii) the galaxy positions are not correlated 
with the shear field, (iii) the galaxy redshifts are not correlated 
with each other or the shear field, and (iv) the galaxies' intrinsic 
ellipticities are uncorrected with each other and the shear, the 
estimators ( [T3) l are unbiased ( [Schneider et al.|2002a| l. 



{i±im) = urn- 



(15) 



The covariance of the cosmic shear estimators determines 
the accuracy to which one can deduce cosmological param- 
eters from these estimators. The covariance c-t-t{§i,§2) = 
Cov(|-t()?i),|'-t(i?2)) of the estimators ([T3| for scales §i and 
§2 can then be split into an ellipticity-noise part c^^l(&i,'&2)^ a 
mixed part c^'^±^(??i, )?2), and a cosmic variance part c^^l('d-i,'0-2): 

c±±(^i,i?2) = cit(t?i,r?2) + cif (t?i,t?2) + c^lli&u&i). (16) 

First, we consider the covariance for a single survey field 
F with area Ap, galaxy density Wg = N„/Ap, and linear dimen- 
sions large compared to -ffi and &2- The ellipticity-noise part c^'^]_ 
vanishes for c+_, and for disjoint bins. Its contribution to the 
auto-variance of the bin with radius {)■ depends on the variance 
crl = (k*'^|^) of the intrinsic ellipticity and the expected effective 
number of galaxy pairs NJd-) - 2n'&A{d-)Apn^ in the bin. 



(17) 



The mixed part c^_^l^ contributes to both the auto- and covari- 
ance of bins: 



m„Ap 



Jo 



(|)5-2«0 -f?!"^,!), 



(18a) 



' see Appendix[c]for a detailed discussion of cosmic shear estimators 
and their noise properties, which also considers a more general estima- 
tor with non-uniform weights for the galaxies. 
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of the cosmic shear covariance reads (see Appendix [C]l 



2crl 

X cos [4^(j?2Mo - '^iM,^,)] , and 



(18b) 



in r% 

/CqAf Jo 

(>?i,i?2)=^rd^i^+(|i?2«()->?i«„l)cos(4^i). (18c) ^ ^nim.m C\Qn ,,Q. (1 
'"'s^^-Jo K^Ap Jo ' Jo 

d03 03f+(03) d 

Jo Jo 



d03 03^+(03) 



Here, denotes a unit vector in the plane with polar angle (p. 

(y) 

Calculating the cosmic variance part c/j! requires knowledge 
about the 4-point statistics of the convergence field. If one as- 
sumes that the convergence field is a homogeneous and isotropic 
normal random field, the cosmic variance parts read ( |Schneider| 
let al.|2002a| l 



4^+(,?i)/7(^2) 



d(^l^+(|6'3Ho -i9-iM^,|) 



d03 03f+(03) d^2^+(|03Mo-1?2M^J) (23) 



+ — ; d6l3 6*3 difx dip 

27rAF Jo Jo Jo 

X^+(|6l3H0 - J?lH^, D^+dflsMO - ^1U^2\) 



4 



Jo Jo 



d<^l ^±()?lM,^,,03«o) 



X I d<^2^±('9-2«:^2'^'3«o) 

Jo 



(19) 



2 + 4^o2^+(03) + ^oV+(03)^+(|03«O - '^IM^, - »?2M^J)1 

with Tjiff) - K^^^+iff) + 1. We are unable to derive equally 



'concise' expressions for the cosmic variance parts c*/i and c 
However, c^il^_ and c"I]_ can be obtained from any given c^ll by an 
integral transformation with the kernel ( fTT) (see Appendix [C|: 

Jr»oo 
M'^§'^G-{^2,»'^)c^^l(M2), (24) 


M\&\G-{^1,&\) I d)?2)?2G-(l?2,'?2) 

Jo 



where 



(^{e,e')^U\e' -e\), 

^_ {e, ff) = - 0\) cos [piO) - <p{0' - 6)] ., 



(20a) 
(20b) 



and Op ~ y/Ap/n denotes the 'radius' of the survey field, a scale 
similar to the field's linear dimensions^ 

In this work, we also consider zero-mean shifted log-normal 
random fields (see Appendix for details) as models for the con- 
vergence field. Such a field k can be obtained from a homo- 
geneous and isotropic normal random field n with mean /i and 
standard deviation cr by the point-wise transformation 



(25) 



Equation ( |23) comprises terms of 2nd to 6th order in If 
only the quadratic terms, which are also present in the expres- 
sion ( [T9| l, and the simplest cubic term in are retained, the cos- 
mic variance contribution reduces to 



4 



X 



Ki0) = exp[n(0)] - Ko, 

which implies a zero-mean shifted log-normal pdf. 



(21) 



Pk(k) 



exp 



[ln(^/^o + 1) + cr^lf 
2o-2 



yJlniK + Ko)cr 



for K > -Ko, 
otherwise. 



(22) 



Choosing the shift ko = exp(yu -i- cr^/2) ensures zero mean. Since 
-a:o marks the lower limit for all possible k, we call kq minimum- 
convergence parameter. 

If the convergence field is well described by a zero-mean ^ 
shifted log-normal random field, the cosmic variance term c^Jl ^^(§) - — ^ f^(i?) 



AQ^Qt, d^i^±()?iH^,,6'3Ho) 

Jo Jo 

r d^2^±('?2«^,,03«o) (26) 
Jo 

+ U^X)U^2)^ Pd03 03^+(03). 

ii()Ap Jo 

In the following, we call the approximation defined by 
Eqs.([T6|-([T9| the normal approximation to the cosmic shear co- 
variance. If Eqs. (|23]l-(|25]l are used instead of Eq. ( [T9| l, we call 
the resulting approximation to the cosmic shear covariance the 
log-normal approximation. If Eq. (|26| is used instead, we call 
the resulting approximation the simplified log-normal approxi- 
mation. 

In general, the convergence field is neither a normal nor a 
zero-mean shifted log-normal field. In case one has a large sam- 
ple of Nr (quasi-)independent realisations of the survey field with 
measured correlations one can estimate the mean using 

the sample mean. 



(27) 



Schneider et al 



1 2002a I and 



Joachimi et al. 



20081 



for simplicity. This is justified for sufficiently large fields, where cor- 
relations on scales larger than the field are negligible, and extending 
the integration boundary of 63 in Eq. l |19[ > to infinity does not affect the 
integral significantly. For small survey fields, however, this may lead 
to an overestimation of the cosmic covariance, which was observed by 
|Sato et aLlj2011b| l in computer experiments with Gaussian convergence 
fields. 



assume 0p w 00 ^jj^j (jjg covariance using the sample covariance, 

1 



(28) 



For many surveys, the survey area is not a single field, but 
consists of several small areas that are far apart from each other 
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on the sky. In this case, the measurements in each field can 
be considered independent from the measurements in the other 
fields. If all survey fields have similar field areas and galaxy den- 
sities, the covariances for the whole survey can be obtained by 
computing the covariances for a single field and dividing the re- 
sult by the number of fields. 

2.3. Estimation of cosmological parameters 

From a set of shear correlation functions ^±()?i), . . . ,^±{&Nhi„,) 
measured at A^bins angular separations, one can form a data vector 



maximum-likelihood estimation with a Gaussian likelihood. In a 
Bayesian analysis with constant priors on the model parameters, 
the likelihood L{7r\d) is proportional to the posterior distribution 
p{n\d). In this case, the above equations also specify how the 
posterior distribution depends on the assumed covariance of the 
cosmic shear correlation functions. For more general prior dis- 
tributions pin), the posterior reads 



Pirn = 



J d^^n L{7r\d)p{7ry 



(36) 



(29) 



In one of its simplest forms, Bayesian model-parameter estima- 
tion from a measured cosmic shear data vector d amounts to a 
maximum likelihood parameter estimation assuming a quadratic 
log-UkeUhood 



In L{n\d) ^--[d- fi{7r)f C^^ [d - ^(;r)] + const. 



(30) 



Here, fiin) denotes the prediction of the cosmic shear correla- 
tion vector for the model characterized by a set of A^p parameters 
;r = (tti , . . . , TTAfp)'. If one uses the normal or log-normal approx- 
imation to the cosmic shear covariance, the inverse covariance 



matrix Cj of the data vector d reads: 



c+ 
c 



with 



= (c..()?/,)?,))|^:;. (31) 



If the cosmic shear covariances C±± are estimated from the 
sample covariances of A^^ realisations, the inverse covariance can 
be estimated by ( |Anderson|2003[|Hartlap et al.|2007 i 



N, - 2Nun. - 2 (c 



(32) 



where the correction factor ensures that this estimate is unbiased 
if the joint distribution of the measured correlation functions is 
well approximated by a multivariate normal distribution. This 
estimate requires A^, > (2A^|,ins + 2). For smaller N^, the sam- 
ple covariance matrix is not invertible, and other methods are 
required to estimate the inverse covariance matrix. 

The parameter set ttml maximizing the (log-)likelihood sat- 
isfies 



= 



dn 



(33) 



A Taylor expansion of In L{n\d) around ttml up to second order 
yields 



\nL(n\d) 
where 



1 



(n - TTuhf C^' ('T - ttml) + const.. 



c ' = 



U,/=i 



1 

d,kl 



dnidiTj 



(34) 



(35) 



differs from the Fisher matrix by additional contributions from 
second derivatives of // when //(^Tml) + d. 

The equations above show how the covariances and con- 
fidence regions of the model parameters depend on the as- 
sumed covariances of the cosmic shear correlation functions in a 



3. Simulations 

We use a numerical approach to assess the performance of pre- 
dictions for the cosmic shear covariance. By ray-tracing through 
the Millennium Run (MR), a large A^-body simulation of cosmic 
structure formation by Springel et al.| ( |200 5 ), we create a suite 
of simulated fields of view with maps of the effective conver- 
gence and shear. From these maps, we measure the convergence 
distribution and the convergence and shear correlations. We then 
estimate the covariance of the shear correlation functions in dif- 
ferent ways: from the measured one- and two-point statistics of 
the convergence using the normal, log-normal, and simplified 
log-normal approximation, and from the sample covariance of 
the shear correlation in the simulated survey fields. 

The MR assumes a flat ACDM cosmology with a mat- 
ter density - 0.25, a baryon density fib = 0.045, and a 
cosmological-constant energy density = 0.75 (in units of 
the critical density), a Hubble constant h - 0.73 (in units of 
100km s"' Mpc"'), a primordial spectral index - \ and a 
normalisation parameter erg - 0.9 for the linear density power 
spectrum. These values also define our fiducial cosmology in the 
later anal ysis. The simul ation employed a customized version of 



GADGET-2 (Springel 2005) with 10'" particles of 8.6 x 10**/?"' M, 



in a cube of 500/z~' Mpc comoving side length, and a comoving 
force softening length of 5/1 ' kpc. 

We employ the multiple-lens-plane ray-tracing algorithm de- 
scribed in [Hilbert et al.| ( [2009| l to simulate gravitational lensing 
observations. The ray-tracing algorithm takes into account the 
gravitational deflection by the dark matter, represented by the 
dark-matter particles of the MR, and the deflection by the stellar 
mass in galaxies (for details, see Hilbert et al.|2008 1 as inferred 
from the galaxy-formation model of De Lucia & Blaizot ( 2007[ l. 

Using the ray-tracing, we generate 64 fields of view. Each 
field has an area of 4 x 4 deg^ (yielding a total area of 1024 deg^) 
and is covered by a regular grid of 4096^ pixels (yielding a res- 
olution of 3.5 arcsec). For each pixel, a ray is traced through the 
MR up to redshift 3.2, and the convergence and shear along the 
ray is recorded. To study the covariance of smaller survey fields, 
we also create 256 fields of 2 x 2 deg^ by splitting the 4x4 deg^ 
fields evenly. 

The convergence and shear information along the rays is 
then used to calculate the effective convergence and shear in 
the simulated fields for a source population with median red- 
shift Zmedian = 1 and a rcdshift distribution ^Brainerd et al.|1996l l 



Pziz) = -3 exp 



-) 



3/2 



, where zq - 



^median 

1.412' 



(37) 



The only source of noise (i.e. uncertainties affecting the esti- 
mation of cosmological parameters) in these ellipticity noise- 
free maps is the cosmic variance (each field represents only a 



5 



Stefan Hilbert et al.: The log-normal approximation to cosmic shear covariances 



measured 
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Fig. 1. The pdf PniK) of the convergence k for sources with me- 
dian redshift Zmedian = 1 ■ The pdf measured from the simulations 
(soUd Une) is compared to a zero-mean normal distribution with 
matching variance (dashed line), and the best-fitting zero-mean 
shifted log-normal distribution (dotted line). 



single finite realization of the underlying cosmology). To cre- 
ate noisy lensing maps that also incorporate the uncertainties in 
cosmic shear surveys arising from the intrinsic shapes of source 
galaxies, we add Gaussian noise with standard deviation crpjx - 

3.5' 



(Tfl -^JlrigApiy^ to each pixel (with pixel area Apix 



^ arcsec^) 



of the simulated shear fields. We assume an intrinsic galaxy el- 
lipticity distribution with standard deviation cr^ = 0.4, and we 
consider surveys with galaxy densities ng — 25 arcmin"^ and 
n„ - lOOarcmin"^. 



4. Results 

4.1. The convergence distribution 

We determine the convergence distribution from the simulated 
fields by binning the ellipticity noise-free convergence values of 
all field pixels into bins of size 0.002. The resulting probability 
density function (pdf) of the convergence distribution is shown 
Fig. [T] For comparison, we also show the pdf of a zero-mean 
normal distribution, whose variance equals the measured vari- 
ance of the convergence, and the pdf (|22]) of the best-fitting zero- 
mean shifted log-normal distribution, whose parameters kq and 
cr where obtained by a simple least-squares fit. Neither the nor- 
mal nor the zero-mean shifted log-normal distribution matches 
the measured convergence distribution perfectly, but (at least vi- 
sually) the zero-mean shifted log-normal distribution fares far 
better than the normal distribution. This finding indicates that 
the log-normal approximation might possibly provide a better 
approximation to the covariance of the cosmic shear correlation 
than the normal approximation. 

While the normal approximation to the covariance of cos- 
mic shear correlations only needs the convergence correlation as 
input, the log-normal approximation also requires the minimum- 
convergence parameter kq. One could, for example, compute the 
convergence for an empty beam to a fixed source redshift and 
use its modulus as kq. Simulations show, however, that there are 
no empty lines of sight in a ACDM universe ( [Taruya et al.|2002] 
IVale & W hite'2003| [Hiibert et alJ2007| l. We thus use the value 
Ko » 0.032 obtained from the fit p2|). 

The minimum-convergence parameter Koiz) for sources at a 
single redshift z, obtained from fits to the measured pdf of the 



Fig. 2. The minimum-convergence parameter /<■() as a function of 
source redshift z (obtained from fits to the pdf of the convergence 
for sources at redshift z). 



Table 1. The parameters kq and cr of the fit ^22) to the con- 
vergence pdf for sources at various redshifts z, compared to the 
measured standard deviation {k^Y^^ of the convergence. 



z 




cr 




0.28 


0.0035 


0.81 


0.0058 


0.51 


0.0100 


0.74 


0.0124 


0.76 


0.0120 


0.70 


0.0199 


0.99 


0.0291 


0.66 


0.0270 


1.50 


0.0542 


0.58 


0.0415 


2.07 


0.0819 


0.53 


0.0547 


3.06 


0.1253 


0.47 


0.0716 


4.18 


0.1608 


0.44 


0.0838 


5.29 


0.1917 


0.41 


0.0922 



convergence k(z), is shown in Fig.|2]as a function of source red- 
shift and listed in Table [T] along with the second fit parameter 
cr and the standard deviation of the convergence. The parame- 
ter cr can be interpreted as a measure of non-Gaussianity of the 
convergence field. The values 0.5 ^ cr ^ 1 indicate a moderate 
degree of non-Gaussianity ( Martin et al.|201 l l. 

The redshift dependence of kq in the range 0.3 < z < 4 is 
well described by 



Ko(z) = O.OOSz + 0.029z^ - 0.0079z^ + 0.00065z^ 



(38) 



The minimum-convergence parameter for sources at redshift z = 
1 is very similar to the value measured from the source redshift 
distribution ( [37] i with median redshift Zmedian - 1- Moreover, 
the minimum-convergence for z = 1 and z = 2 are in good 
agreement with the corresponding values found by Taruya et al. 
(j^02|. 

4.2. The convergence and shear correlations 

From each ellipticity noise-free simulated 4x4 deg^ field, we 
estimate the convergence correlation ^^(6) and the cosmic shear 
correlations ^-t using Fast Fourier Transform (FFT) techniques 
employing zero-padding and book-keeping of the number of 
contributing pixel pairs to properly account for the non-periodic 
field boundaries (see Appendix[D]for more details). The correla- 
tion estimates from all fields are then used to compute the mean 
correlations and their uncertainties from the field-to-field vari- 
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Fig. 3. The convergence correlation ^k{^) (squares) and the cos- 
mic shear correlations (diamonds) and ^-{ff) (triangles) as 
functions of the angular scale i? for sources with median redshift 
^median = 1 measured from the simulations. Error bars indicate 
uncertainties on the mean estimated from the field-to-field vari- 
ance of the simulated fields. Also shown are the analytic fits (see 
text) to the correlations ^+()?) (dashed line) and ^-(fl) (dotted 
line). 



Or simulations show furthermore that the normal approxi- 
mation completely fails to reproduce the measured c'll and cIl 
on all scales considered. As Eqs. (flOll, ( |25| l, and (|24| show, the 
shear correlation f_ and the cosmic variance parts c_l_ and c__ 
at given scales depend on the statistical properties of the conver- 
gence field on all smaller scales. These scales, of course, include 
the non-Gaussian scales < 10 arcmin, where the normal assump- 
tion fails. 

In contrast to the normal approximation, the log-normal ap- 
proximation yields estimates of the cosmic variance terms c^^ 
that are very similar in shape and magnitude to the measured cos- 
mic variance contribution. Differences are noticeable on scales 
below a few arcmin, where the log-normal approximation over- 
estimates the covariance. Furthermore, the log-normal approx- 
imation underestimates the covariances involving ^_ on scales 
> 10 arcmin. The same holds for the covariances obtained from 
the simplified log-normal approximation (not shown), which are 
almost identical to the covariances from the log-normal approx- 
imation. 



4.4. Comparison to empirical fits for tine cosmic shear 
covariance 



The resulting convergence and cosmic shear correlations are 
shown in Fig. [3] As expected, the measured shear correlation 
f+ is almost identical to the measured convergence correlation 

Numerical tests verify furthermore that the measured corre- 
lations and ^_ satisfy the relation ([TOf well within the error 
bars. 

Computation of the covariance of the cosmic shear correla- 
tions in the (log-)normal approximation involves integrals over 
the shear correlations. To facilitate the integration, we approxi- 
mate the measured correlation functions by analytic expres- 
sions 



1 



and 



(1 + §lT^f 



(39a) 



(39b) 



with parameters ai, 02, b^, b2, b^, tj, T2, and T3 obtained from a 
weighted least-squared fit. The resulting expressions fit the mea- 
sured correlations very well, as Fig.|3]shows. 

4.3. The covariance of the cosmic shear correiation functions 
We estimate the cosmic variance contribution cTl to the cosmic 



shear covariance from the sample covariance ( [27| i of the mea- 
sured shear correlation functions |^(i?) in the ellipticity noise- 
free simulated fields. Furthermore, we calculate estimates for 
Cj_" within the normal approximation and the log-normal ap- 
proximation, using the fits p9a| i as inputs for the shear correla- 
tion functions. The estimates from the sample covariance of the 
2 X 2deg^ fields are compared to the corresponding estimates 
based on the normal and on the log-normal approximation in 

Fig-El 

As is already known from earlier studies (e.g. |Semboloni| 
et al.|2007 Sato et al. 201 Vo) , the normal approximation grossly 



underestimates the covariance on scales < 10 arcmin. The 
reason for this failure is that the convergence and shear fields are 
very 'non-Gaussian' on such small scales, where the growth of 
matter structures is very non-linear 



The failure of the normal approximation on small scales has mo- 
tivated Semboloni et al. (200T) and Sato et al. ( 201 lb ) to develop 
empirical corrections. Figure |5] shows the covariances resulting 
from applying the correction factors they proposed for sources 
at redshift z = 1 to the cosmic varianc e part c^ j in the normal 



approximationplThe estimate based on Semboloni et al. (2007i 

Tr) 11 u.. „ , f„„._ .1.' 



overestimates c 'l on small scales by a large factor, whereas the 
corrections according to |Sato et al.| ( |201 Fb) ! yield estimates that 
underestimate cill by 20-50%. 



Semboloni et al. (2007i and Sato et al. (201 lb I do not pro 



vide corrections for the cosmic variance parts or c'lL We 
thus employ Eqs. ( |25] l and ( [24]l to compute c'/i and c"Il from 
the predictions for As Fig.pjshows in comparison to Fig. [4] 
the estimates based on Semboloni et al. ( 2007| l fail to reproduce 
the estimates from the lensing simulations both quantitatively 
and qualitatively. On most scales, the estimate is too large by 
orders of magnitude. Around the transition scale between the 
Gaussian and non-Gaussian regime at 10 arcmin, the estimated 
covariances show artificially sharp drops. 



The estimates based on 



Sato et al. (201 lb I perform better 



However, they still overpredict the covariances on small scales 
by a factor of a few, and underpredict the covariances on large 
scales. Since the empirical corrections were designed to yield 
good estimates of c'll for the range of scales and cosmologies 
considered by |Semboloni et al. ( 2007 [ ) and [Sato et al. (2011b), 
this indicates that we have used the corrections outside their 
range of applicability. For the cosmology and scales considered 
here, the log-normal approximation appears to perform better 
than the empirical corrections. 



The corrections by Sato et al. 1 2011b 1 assume explicitly that the 
normal prediction Eq. \\9\ is computed with ^ = 00 (as in |Schneider| 
et al. 2002a). We thus apply their correction factor to the normal predic- 
tion computed assuming % = 00 instead of 9-g = y/ApJn. Moreover we 
extend the range of the fit correction F defined in eq. (Bl) to all scales 
where F > 1 to avoid discontinuities in the covariance. 
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Fig. 4. The cosmic variance contributions c^l({)-i , i?2) to the cosmic shear covariance as functions of angular separations i?i and )?2- 
Shown are the contributions c^^l (left column), c^^l (middle column), and c^i''! (right column) for 2x2 deg^ fields measured from 
the lensing simulations (top row), computed employing the normal approximation (2nd row), and the log-normal approximation 
(bottom row). 



4.5. Positive-semidefiniteness of the cosmic variance 
contribution 

In the absence of other noise (e.g. ellipticity noise and other 
sources not considered here), the cosmic variance part of the 
cosmic shear covariance represents the full data covariance. 
This implies that the actual cosmic variance part of the pa- 
rameter covariance matrix as well as its components c^:^l and 
c__ are positive-semidefinite regardless of the statistical prop- 



erties of the convergence field. Positive-semidefiniteness of the 
data covariance is, moreover, an essential requirement for a 
sound Bayesian parameter inference employing a quadratic log- 
likelihood. Thus, a highly desirable property of any approxi- 
mation to the cosmic variance contribution is that it yields a 
positive-semidefinite cosmic variance part, so that even in cases 
of high galaxy densities and low ellipticity noise, the full data 
covariance remains positive-semidefinite. 
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The cosmic variance contributions c"ll{-d-\,§2) (left column), cfl{-d-\,§2) (middle c olumn), and c^^l(i9i,§2) (right column) as 



Fig. 5. 

functions of angular separations )?i and §2 for 2 x 2deg^ fields computed from the fit by Semboloni et al. (2007 top row) and the 
fit by|Sato et aL](|2011b[ bottom row). 



In the cases studied here, we encounter a serious problem of 
the corrections suggested by Semboloni et al. ( 2007) 1 and |Sato| 
et al. (201 Ibj): Both fail to yield positive-semidefinite covari- 



ance parts cfl and dlL As a consequence, the cosmic variance 
part of the data covariance matrix is not positive-semidefinite, 
but has significantly negative eigenvalues. This may lead to a 
non-positive data covariance matrix Cd, a non-positive parame- 
ter covariance matrix C^j, and to a breakdown of the parameter 
estimation. 

Apparently, the empirical corrections do not contain any 
mechanism ensuring positive-semidefiniteness of the data co- 
variance. Indeed, it seems difficult to devise a valid covariance 
matrix of correlation functions without strong guidance from, 
e.g., analytic models. 

In contrast, both the normal and the log-normal approxima- 
tion 'should' yield positive-semidefinite covariance matrices by 
construction, since the resulting matrices are covariance matri- 
ces of correlation functions of a (log-)normal field Within the 
simplified log-normal approximation, the cosmic variance part 
of the covariance matrix is a sum of the cosmic variance part 
resulting from the normal approximation and another positive- 
semidefinite matrix, and thus is also expected to be positive- 
semidefinite. However, due to the employed approximations for 



* This certainly holds if the cosmic shear correlations assumed in 
the computation of the covariances are valid correlation functions of a 
(log-)normal field, which we checked for our calculations. 



finite fields of view and numerical inaccuracies, the normal and 
(simplified) log-normal approximations may also yield indefi- 
nite cosmic variance parts with negative eigenvalues. These neg- 
ative eigenvalues are found to be essentially consistent with zero, 
indicating the presence of linear constraints on the data vector 
[e.g. stemming from the relations (|9]) and {T0\], and of much 
smaller magnitude than the negative eigenvalues encountered 
when using the empirical fits. Setting the negative eigenvalues 
to zero 'by hand' does not change the inferred posterior distri- 
butions of the parameters significantly for the case of the normal 
and log-normal approximation, as shown in the following sec- 
tion. 



4.6. Error estimates for cosmological parameters from 
cosmic siiear 

To study how the approximations to the cosmic shear covariance 
affect the inferred accuracy of cosmological-parameter estima- 
tion, we consider the scenario that cosmic shear data is used to 
constrain the parameters and erg of a flat ACDM cosmol- 
ogy. To keep the discussion simple, we choose priors that are 
constant for Q.^ e [0.15,0.35] and erg € [0.7, 1.1] and zero oth- 
erwise. Furthermore, we assume the other cosmological param- 
eters known to take the fiducial values: h - 0.73, Q.\, - 0.045, 
and rts - I. 

We assume the cosmic shear correlation functions ^+()?) and 
^_()?) are measured in 20 logarithmically spaced bins in the 
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Fig. 6. Joint posterior distributions of the mean matter density Qm and the matter power spectrum normalization erg of a flat ACDM 
universe inferred from a cosmic shear survey with Wg - 25 arcmin"^ (assuming flat priors for and erg, and all other cosmological 
parameters known). Shown are the 68% (inner) and 95% (outer) confidence contours for a small survey in 5 fields of 2 x 2 deg^ (left 
panels), a survey in 6 fields of 4 x 4deg^ (middle panels), and a large survey in a 120 x 120 deg^ field (right panels) obtained when 
using the covariance matrix measured from our lensing simulations (solid lines), the normal approximation (top, dashed fines), the 



log-normal approximation (top, dotted lines), the empirical corrections to the normal approximation by Semboloni et al. (2007 
bottom, long dashed lines), and the corrections by |Sato et al.| ( |20lTb] bottom, long dash-dotted lines). 



range 0.2arcmin < § < 120arcmin. The model predictions 
//(Qm, CTg) for the shear correlations are computed using nicaea 
( Kilbinger et al.|[2010| l. The data vector d containing the mea- 
sured shear correlations is assumed to coincide with the pre- 
dicted values for the fiducial cosmology with Q.^ - 0.25 and 
(Tg = 0.9. 

As first example, we study the parameter constraints ob- 
tained from a small survey with n„ - 25 arcmin"^ in 5 fields 
of 2 x 2 deg^ (i.e. a total area of 20 deg^) similar to the Deep 
Lens Survey (see, e.g., Kubo et ar||2009 1. We estimate the in- 
verse covariance matrix directly from the sample covariance of 
the simulated noisy shear maps using Eq. ( [32] l, and combine the 
estimate with the nicaea predictions to compute the likelihood 
and the posterior. 

In the left panels of Fig.|6] the resulting confidence contours 
in the (iln,, cr8)-plane are compared to the contours obtained 
from the various approximations to the cosmic shear covariance. 
The normal approximation substantially underestimates the size 
of the confidence regions, in particular in the direction perpen- 
dicular to the major degeneracy. In contrast, the contours based 
on the log-normal approximation are very similar to the con- 
tours based on the measured covariances. The same holds for 
the simplified log-normal approximation (whose contours are 



not shown, since they are almost identical to those from the log- 
normal approximation). Only a slight tilt of the contours based 
on the (simplified) log-normal approximation against the con- 
tours based on the ray-tracing is visible. The corrections pro- 
posed by Semboloni et al.| (j2007) yield confidence regions that 
are much larger than the regions infeiTed from the simulations, 
and are strongly tilted. The confidence regions based on the cor- 
rections proposed by |Sato et aL] ( |201 lb| l are noticeably smaller 
than the regions inferred from the simulations. 

As second example, we consider a survey with rig = 
25 arcmin"^ in 6 fields of 4 x 4 deg^ (i.e. a total area of 96 deg^) 
similar to the Canada-France-Hawaii Telescope Wide Synoptic 
Legacy Survey (Hoekstra et al.^2006) . The resulting confidence 
contours are shown in the middle panels of Fig.|6] As in the first 
case, the normal approximation underestimates the size of the 
confidence regions, whereas the (simplified) log-normal approx- 
imation yields confidence regions similar to those infeiTed from 
our lensing simulations. 

As can be seen in the right panels of Fig. |6] the normal ap- 
proximation substantially underestimates the size of the con- 
fidence regions also for very large surveys like the planned 
Large Synoptic Survey Tel escope (LSST) survey (LSST Science 
Collaborations et al.|2009 1 or the Euclid survey (Refregier et al. 
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|2010| l. Even for such large surveys, much cosmological informa- 
tion is contained in the shear correlations on small scales, which 
are not well described by the normal approximation. The (sim- 
plified) log-normal approximation, describing the small scales 
better, yields confidence regions in much better agreement with 
the regions inferred from our lensing simulations (using a scaled 
version of the cosmic covariance measured in the 4x4 deg^ fields 
in combination with analytic expressions for the mixed and el- 
lipticity noise). The contours based on Semboloni et al. ( 2007| l 
appear much too large and strongly tilted, whereas the contours 
based on| Sato et al.| ( 20lTbl ) appear too small. 

In the cases just discussed, both the cosmic variance and the 
intrinsic ellipticity noise are important sources of error for the 
inferred parameters. To investigate the case that the cosmic vari- 
ance is the dominant source of parameter uncertainty, we con- 
sider again a survey in 5 fields of 2 x 2 deg^, but with a higher 
galaxy density rig - 70 arcmin"^. The left panels of Fig. [vjcom- 
pare confidence contours estimated from the simulations and 
the contours obtained from the approximations. The normal ap- 
proximation underestimates the credible parameter regions even 
more than in the afore discussed cases. The (simplified) log- 
normal approximation still yields confidence regions compara- 
ble in size to those estimated from the simulations, but the shape 
difference is more apparent (this might become problematic in 
cases where the compatibility of different cosmological exper- 
iments is judged by the overlap of their parameter confidence 



regions). The confidence regions based on the fit by |Semboloni 
etaL| ( [2007l l are of similar size, but strongly tilted. 



When using the fits proposed by Sato et al. ( 201 lh\ , we do 
not obtain any reasonable confidence contours for the posterior 
Both data and parameter covariance matrices are indefinite, and 
the likelihood features a saddle point at the fiducial parameter 
values instead of a maximum. The problem can be traced back 
to the cosmic variance part of the parameter covariance matrix, 
which is indefinite. 

One may approach the problem of negative eigenvalues in 
the parameter covariance matrix stemming from an indefinite 
cosmic variance part as follows: Using its eigensystem decom- 
position, the cosmic variance part of the covariance matrix can 
be specified by an orthogonal set of eigenvectors, which describe 
the principal directions of scatter in the data vector due to cosmic 
variance, and the corresponding eigenvalues, which quantify the 
scatter in these directions. Directions with vanishing scatter in- 
dicate linear constraints, which confine all possible data vectors 
to a lower-dimensional subspace of the full data vector space (in 
the absence of other sources of scatter). Directions with negative 
eigenvalues, i.e. 'negative scatter', do not make sense. One may 
presume, however, these negative eigenvalues stem from numer- 
ical inaccuracies in the employed method for computing the cos- 
mic variance, and rather indicate directions with very small pos- 
itive or vanishing scatter 

To ensure a positive-semidefinite data covariance matrix 
from the various approximations, we modify the cosmic variance 
part by replacing any of its negative eigenvalues in its eigensys- 
tem decomposition by zero. The resulting confidence contours 
are shown in the right panels of Fig. |7] There is no visible dif- 
ference in the contours between the left and right panel for the 
normal and log-normal approximation, even though their cos- 
mic variance parts have a few negative eigenvalues (which are 
of tiny magnitude). In contrast, the confidence contours based 
on the fit by Sembol oni et al.| ( [2007] ) change drastically if the 
data covariance matrix is modified in the described way. Only 
after applying the modification, we obtain confidence contours 
for the correction suggested by |Sato et al.| ( |201 lb| l. The param- 



eter confidence regions computed from the modified empirical 
corrections are much smaller than the confidence regions esti- 
mated from the simulations. 



5. Summary and discussion 

Accurate estimates for the covariance of cosmic shear correla- 
tion functions are essential for reliable estimates of the errors 
on the cosmological parameters inferred from cosmic shear sur- 
veys. In this work, we developed two approximations to the 
cosmic shear covariance based on the statistics of log-normal 
random fields. We used numerical simulations of cosmic shear 
surveys to assess the performance of this log-normal and sim- 
plified log-normal approximation and to compare them to the 
widely used normal approximation to the cosmic shear covari- 
ance ( [Schneider et al. 2002a ). 

We find that the normal approximation to the cosmic shear 
covariance substantially underestimates the inferred parameter 
confidence regions, in particular for surveys with small fields 
of view and large galaxy densities, but also for very large 'all 
sky' surveys (like the proposed Euclid or LSST lensing surveys). 
The log-normal approximation yields much more realistic con- 
fidence regions at the price of slightly more complicated expres- 
sions for the cosmic variance contribution to the cosmic shear 
covariance. In contrast, the simplified log-normal approximation 
is as simple as the normal approximation, yet appears as accu- 
rate as the log-normal approximation. Moreover, the simplified 
log-normal approximation is simpler than several proposed ap- 
proximations based on halo models (e.g. Takada & Jain 2009} 
Piel orz et al.j2 010 ). The simplified log-normal approximation is 
also more general than approximations based on empirical fits to 
numerical simulations ( [Semboloni et al.|2007[|Sato et al.|201 lb| l. 

A particular advantage of the normal and (simplified) log- 
normal approximation over the empirical fits suggested by 
Semboloni et al. (2007) andlSato et al. (201 lb) is that the former 



yield positive-semidefinite data covariance matrices by construc- 
tion, whereas the latter may fail to do so. A positive-semidefinite 
data covariance matrix is, however, essential for a sound param- 
eter estimation employing a quadratic log-likelihood. 

A disadvantage of the (simplified) log-normal approxima- 
tion in comparison to the normal approximation is the need 
for providing a minimum-convergence parameter kq. The value 
of this parameter depends on the assumed cosmology as well 
as the source redshift distribution of the weak-lensing sur- 
vey. However, computing /<■() does not require more effort than 
computing the expected shear correlation functions (which are 
needed in any case). Estimates for /to from simulations, for ex- 
ample, require much fewer realisations than estimates for the full 
cosmic shear covariance. 

Because of its comparable simplicity and much better ac- 
curacy, one should consider the simplified log-normal approxi- 
mation in favour of the normal approximation, e.g. for the pa- 
rameter estimation from current surveys, or for parameter error 
forecasts for future surveys. For the analysis of observed cosmic 
shear data from future large and deep surveys, however, better 
descriptions of the cosmic shear covariance will be required. 

In this work, we concentrated on the covariance of the cos- 
mic shear correlations ^-t and the resulting errors on cosmolog- 
ical parameters. In future work, one should adapt the (simpli- 
fied) log-normal approximation to derive covariances for other 
lensing two-point statistics (e.g the COSEBIs introduced by 
Schneider et al. 2 010). The (simplified) log-normal approxima- 
tion could also be generalized to provide covariances for tomo- 
graphic shear surveys. Furthermore, one could consider the ap- 
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Fig. 7. 68% (inner) and 95% (outer) confidence contours of the joint posterior distributions of the mean matter density and the 
matter power spectrum normaUzation erg of a flat ACDM universe inferred from a cosmic shear survey with «g - 70 arcmin"^ in 5 
fields of 2 X 2 deg^, obtained using the covariance from our lensing simulations (solid lines), the normal approximation (top, dashed 



lines), the log-normal approximation (top, dotted lines), the empirical corrections to the normal approximation by Semboloni et al. 
( |200 7' bottom, long dashed lines), and the corrections by |Sato et al.| ( |2011b ' bottom, long dash-dotted lines). For the left panels, 
the cosmic covariance parts of the data covariance matrices have be used as given by the approximations. For the right panels, the 
cosmic variance parts have been modified where needed to ensure positive-semidefiniteness (see text for details). 



proximation of a log-normal convergence field for predictions of 
higher-order shear correlations and their covariances. 

One should also investigate cosmic shear covariances for 
convergence fields that are more general transformations of 
normal random random fields (as considered, e.g., by Das & 
|Ostriker|2006HJoachimi et al.|201 l| l. Such work should also take 
into account the non-Gaussianity and cosmology-dependence of 
the cosmic shear likelihood (Eiflei^ ^ al.||2009 ; Hartlap et al 



2009] [Schneider & Hartlap„2009, ,Sato et al.||2011a, ,Joachimi 
& Taylor||201 1 1. Finally, one should investigate further to what 
extent the matter density field or the convergence field can be de- 
sciibed by a transformed normal random field (see, e.g, Neyrinck 
[eFal . 2009; Seo et al. 2 0TT]|Yu et al.|2011HJoachimi et al.|20lT] l 
and what the physical reasons behind the successes and limits of 
such a description are. 
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Appendix A: Normal random fields 

Consider a non-degenerate homogeneous and isotropic normal random field n in D dimensionsj^This implies that the probability 
density function (pdf) of the field values n - nix) at any single point x e MP is given by 



p„(n\n,a-) 



V27ro- 



exp 



2fr2 



(A.l) 



where ^ and cr^ denote the position-independent mean and variance of the random variable n{x). Moreover, the joint pdf of the field 
values M = («i , . . . , Wiv)' = {n{xi), . . . , n{xN)) at a set of mutually distinct points Xi, . . . ,Xn is given by the pdf of a non-degenerate 
multivariate normal distribution: 



p„in;fi,r) 



1 



(2;r)^/2 VdetZ 



exp 



1 



(n-fi)'Z-\n-fi) 



(A.2) 



Here, fi - {n) - jjl with 1 = (1, . . . , 1)' in accordance with the requirement of homogeneity. The covariance matrix Z is a positive- 
definite real symmetric matrix, whose elements Z,j are given by the two-point correlations ^„jj of the field n. Since the field is 
homogeneous and isotropic, the two-point correlation ^„ ,j of field values at points jc, and Xj depends only on their separation, i.e. 



^,7 = ^",'7 = ^n(\Xi - Xj\), 



(A.3) 



where denotes the correlation function of the field n at separation §. In particular, the diagonal elements of Z are given by the 
variance cr^, i.e. Zu = ^(0) = cr^. 

Using the joint pdf ( |A.2| i, the expectation value {f(n)} of a function /(«) depending on the field values n - («(jci), . . . , n(xN))'' 
at positions x-i, . . . ,Xn can be written as 



</(«))= d^n 



R« (2;r)W2 VdetZ 



exp 



--(n-fifr^n-fi) 



fin). 



(A.4) 



The moment-generating function for the field values « = (ni , . . . , n^)' at mutually distinct points Xi, . . . ,Xn is defined as the 



map 



M„ : ^ R : < M M„{t) = (exp(i'«)) . 
For a normal random field, 

M„{t) = /exp(^'«)\ = f d^M ^ exp 

^ ' Jr« (27r)'v/2 VdetZ 

Jr' 



(A.5) 



--(n-nyi '(n-n) + t'n 



(2n)^/^- VdetZ 



— (n-fi- Zf)' Z ' (n-fi-Zt) + t^n + -t^Zt 







= exp 


_ 2 



(A.6) 



Using M„, the nii ■ ■ ■ m^-point correlation function composed of powers nii, 
by 



For example. 



d" 



"Mnit) 



Qf" ■ ■ ■ dt"" 



t=0 



Qiny-mM 



t^u + -f'Z/ 
^ 2 



, m^v of the field values n i , . . . , can be computed 

(A.7) 



{11x112) = yU + ^„,n, 
{111112113} = jJ^ +H {^n,n + ^«,13 + ^n,2i) , 
{111112112114} - fl" + {^nM + ^n,13 + ^n,14 + ^n,23 + ^n,24 + ^iim) + ^n,12^n,34 + ^n.I3^n,24 + ^n,14^n,23- 



(A.8) 
(A.9) 
(A. 10) 



In a similar manner, one can compute expectation values for products of powers k\,...,kN of exponentiated field values 
exp(ni), . . . ,exp(nAr): 



(exp(ni)'" ■ ■ ■ exp(niv/"^ = ^exp(A:'«)^ - exp 

_ t'(A:-l)/2 pf kikj 



r ^ / 2\ 

= exp ^ ^ -H y -H ^ kikj^nilXi - Xj\) 

i=l ^ ' i<j 



(A. 11) 



where A - exp (fj. + 0-^/2^ k - {k\, . . . , k^y, 770 - exp (cr^j, and 77,^ = rjQxi - x j\) with //(r) = exp [^„(r)]. In particular, 

(exp(ni)) = ^ = exp (ju -i- 0-2/2) . (A. 12) 

^ Readers interested in a rigorous definition of a normal random field should consult the standard fiterature on measure theory. 
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Appendix B: Zero-mean shifted log-normal random fields 

A zero-mean shifted log-normal homogeneous and isotropic random field z can be defined via 

z : ^ R : jc z{x) = exp [n(jc)] - A, (B. 1) 

where n is a homogeneous and isotropic normal field with mean and variance cr^, and the shift A - exp(^yu + to ensure 

vanishing mean for z- The one-point distribution of the random field z reads 



Pziz; A, cr) = 



1 r [\n{zlA+\) + (T^I2f\ 

exp< T > for z > -A, and 



^l2^{z + A)a I 2(^2 j ~ ' (B.2) 

otherwise. 



The negative of the shift A marks the lower limit for possible values for z- 

The A^-point correlation of z can be deduced from the two-point correlation of the underlying normal field n as follows. First 
note that an mj ■ ■ ■ w^-point correlation function of z composed of powers nii,. . . ,mt^ of the field values z\, - ■ ■ ,Zn 31 N mutually 
distinct points Xi, . . . ,Xn can be written as: 



where k = (ki, . . . , k^f and m = (nii, . . . , nif^y-. Combining the results of the Eqs. ( |B.3[ ) and ( |A.ll| i, one obtains 

mi / , niN I \ 

{zT ■ ■ ■ = ^-'^ (;■) ■ ■ ■ 2 (';;:)(-i)*'""^"'?o n ^^-^^ 

For example, 

<ZiZ2> = ^'[7712 - 1] = i'{exp [U\xi - xj\)] - 1} = U\Xi - Xj\), (B.5) 

{Z1Z2Z3} = /I^['7i2?7i3'723 - '712 - rjn - '723 + 2], and (B.6) 

<ZiZ2Z3Z4) = A\rinrinJiumimAVu - miVnm-i - 1nJ]Hm4 - '7i3?7i4'734 - m^m^lM + '7i2 + ?7i3 + ^714 + ?723 + ?724 + ^734 - 3]. (B.7) 



From Eq. ( |B.5| l follows that the correlation functions ^^(r) and ^„(r) are related via 
^,(r) = ^Vp[a'-)]-l} « ^„(r) = lnU-2^,(r)+l]. (B.8) 
The correlation function ^-(r) and the function //(r) = exp are related via 

^,(r) ^ A\i](r) - I] « 77(r) = i-2^,(r) + l. (B.9) 
Using these relations, one obtains, e.g. 

(zizi) = ^.,12, (B.IO) 

{Z1Z2Z3) = fe,12^z,13 + ^z,12^z,23 + ^z,13^z,23] + ^"^^z,12fz,I3fz,23, and (B. 11) 

{Z1Z2Z3Z4} - ^zM^zM + ^z,13^z,24 + ^z,14^z,23 + A'^^ [^,,i2^z,13^z,14 + ^z,12^z,13ft,24 + ^z,I2^z,I3^z,34 + . . .] 

-4 -6 -8 (B.12) 

+ ^ [^z,12^z,13^z,14^z,23 + ■■■]+ A [^jj2^z,13fz,14^z,23^z,24 + ■■■]+ A ^2j2^z,I3^z,14^z,23^z,24^z,34- 



Appendix C: Statistical properties of the shear-correlation estimator 

In this Section, we give a detailed discussion of approximations to the noise properties of the commonly used shear-correlation 
estimators ( [T3| ). The discussion loosely follows the real-space approach of [Schneider et al.] ( |2002a| l. In contrast to previous works, 
we do not restrict the discussion to Gaussian shear fields. Instead, we first derive expressions for the covariance of the shear- 
correlation estimators for the case of a general convergence field, and then discuss the special cases of normal and zero-mean 
shifted log-normal convergence fields. 
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C.1. Basic assumptions and definitions 

To obtain analytical expressions for the statistical properties of the shear-correlation estimator, we consider an ensemble of weak- 
lensing surveys with A^g galaxies in a rectangular field F with area Ap, which implies a galaxy number density rig - Ng/Ap. To keep 
the computations simple, we make the following additional assumptions: 

- The number of galaxies is large, i.e. A^g » 1. 

- Periodic boundary conditions are assumed for the distortion field (i.e. we neglect effects of field boundaries and distortions 
caused by structures outside the survey area). 

- The shear field y is generated by a convergence field k in the survey field (i.e. we neglect B-mode shear). 

- The convergence field Kisa realisation of a statistically homogeneous and isotropic random field (we ignore that the assumption 
of statistical isotropy is in conflict with the assumption of periodic boundary conditions). 

- The galaxy positions . . . ,0*^* are uniformly and independently distributed in the survey field (i.e. we neglect correlations 
between the galaxy positions). 

- The galaxy redshifts z*'', . . . ,2*^^' are identically and independently distributed according to a redshift distribution p-(z) (i.e. we 
neglect correlations between the galaxy redshifts). 

- The observed image elliptic! ty y''' = e^'' + y*'' of each galaxy is a sum of its intrinsic ellipticity e^'' and the shear y^'' - 'y(6^'\z^'^) 
at the galaxy's position 0*'' and redshift z''*. 

- The intrinsic ellipticities e*'\ . . . , e*^^ are identically and independently distributed and are not correlated with the shear field 
(i.e. we neglect intrinsic ahgnments). The joint distribution of the two components ej'^ and of each intrinsic ellipticity e*'* is 

isotropic with vanishing mean and variance cr^ /V2 per dimension (i.e. we neglect complications due to non-isotropic intrinsic 
ellipticity distributions caused, e.g, by insufficient PSF corrections). 

The Fourier transform f{{) of a real function f(6) on the survey field is defined by the relations: 

/(f) = ^ f(f0e-'''>m « m^Yt^'-'m. (cd 

Expectation values (/) of observables f(0^^\ . . .) are obtained by computing the following (formal) average: 

<./> = (f\z,.,. = mf)el)e). (C.2) 

Here, (f)^ denotes the ensemble average over the reaUsations of the convergence field, {/)^ denotes the average over the intrinsic 
ellipticities, 

(f\ = J iz^'^P^iz''^) ■ j dz('^.);,,(z(^^))/(0('), . . .) (C.3) 
denotes the average over the galaxy redshifts, and 

{f\ = ^ r i^^c) . . . r . . .) (C.4) 

Ap° -"F Jf 

denotes the ensemble average over the galaxy positions. We assume that the order of the averages does not matter for all observables 
of interest. 

The tangential component jriS, if , z) and cross componen t jx{0, z) of the shear y(0, z) at position and redshift z with respect 
to the direction are defined by (e.g. |Schneider et al.|2002b| 

yA0, z) = -51 {y(0, z)e-2'*"W) and y^iO, z) = -3 {yi0, z)^-^''''^'*^) , (C.5) 

where (f{&) denotes the polar angle of the vector §. 

We denote the tangential and cross component of the shear y*'^ = y(0*'\z*'') at position 0^'^ and redshift z*'^ of a galaxy / with 
respect to the direction = (f(0^^'' - 0*'*) towards a galaxy j at position 0^^^'' by 

y<;'^' = y^(6/», 6l(^) - 0^'\ z('>) and y'^i^^ = yx(6l<'\ 0^^^ - 6/«, z^''). (C.6) 
Similar definitions are made for the tangential and cross component of the intrinsic and observed image elhpticity: 

= -% (e('>e-2V'-'')) , = _3 (^(0e-2V-') , (C.7) 

y'i'^'^ = -% (y(''e-2'^"'^') = y^»^ + and y*;-^' = -3 = y^'^> + e^'^l (C.8) 

Furthermore, we define the source-redshift averaged convergence k(0) and shear y(0) as in Eqs. ([3]) and (|4]). From the assumption 
that the shear field is generated by the convergence field follows that the convergence and the Cartesian shear components are related 
in Fourier space by: 

yi(e) ^ cos[2ipit)]K(e) and y2it) ^ sm[2(pit)]K(f). (C.9) 
We define the tangential and cross components of the effective shear as in Eqs. (|5]l and (|6]l. With these definitions, 
yAO, &)yr(0 + M + yxi0, ^)7x(0 + n(O)yi(0 + #) + yi{0)y2iO + and (C. 10) 

yA0,^)yA0 + ^,^)- yAe,&)yA0 + i?,i?) = cos[4^(#)][yi(0)yi(0 + ,?) - y2(0)y2i0 + &)] 

+ sin[V(#)][yi(0)y2(0 + #) + y2(0)yi(0 + i?)]. 



16 



(C.13) 



Stefan Hilbert et al.: The log-normal approximation to cosmic shear covariances 
C.2. Correlation functions and tlieir relations 

There exist well-known relations between the two-point correlations of the convergence and the shear, e.g. the Eqs. (|9| and ( [TO] l. 
We wish to derive similar relations between the four-point congelations of the convergence and the four-point correlations of the 
shear, which appear in the covariance of the cosmic shear estimators ( [T3] l. Achieving this task for general convergence fields 
requires relations between the convergence and the shear that hold not only for the ensemble mean, but for every realisation of the 
convergence field. 

For every realisation, we define 'empirical' coiTelations based on spatial averages: 

UO) = r i^O'KiO'MO' + 0) and UO) = f [TtC^', e)y^{e' + 0,0)± y^{0\ 0)y^ie' + 0, 0)] . (C.12) 

'^F Jf Jf 

From the convolution theorem and the relations ( |C.9| l, ( |C.10| l, and ( |C.11[ ) follows that for every reaUsation, 

UO) = -^ r i^O' [yi(0')7i{0' +0) + J2{0')y2{e' + 0)] = V [co^[2ip{t)f + sin[2^(f)]'} r(t)K{t) = V e'"r(fM^) 
Af Jf V V 

= UO), and 

= ^ r l'^'{cos[4^(e)] [y,(0)yi{0' +0)- y2{&)yi(0' + 0)] + sm[M0)] [yi(0')72(O' + 0) - yi(&)yi(0' + 0)]} 

'^F Jf 

= " {cos[4^o(6')] cos[2^(f)]^ - cos[4,^(6')] sin[2^(^)]^ + 2 sin[4,^(0)] sin[2^o(^)] cos[2^(/')]} ~K{t)K{t) 

= y e'''" co5[Aip{{) - Aip{0)\K{t)K{{) = y e'''-^ cos[4<^(f) - Aip{0)\ — f i^ffe-'^"' ^A0')- 
This can be written as 

UO)^ {a^0'G^{0,0')UO') (C.15) 

Jf 

with the kernels 

GA0,0') = 6o(0-0') and G40,0') = — V e'^'^"''^' cos[4-ip({) - Aip(0)], (C.16) 

Af V 

where 6d denotes the (in this case two-dimensional) Dirac delta 'function'. The relations ( |C.15| l between the empirical convergence 
and shear correlations are valid for every realisation. They will provide us with the tools to translate the four-point correlations of 
the shear into four-point correlations of the convergence. Of course, they can also be used to derive the familiar relations between 
the two-point correlations of the shear and the convergence. 

The ensemble-average convergence and shear correlations are defined by (e.g. Blandford et al.|1991[|Kaiser|1992 l: 

Um = ilm = {<0M0 + i?)) and = {^^m = {yA0, ^)y^{0 + + {yA0, ^)yx(0 + ^, ^)) . (C.17) 

The statistical isotropy of the convergence field implies that its correlation and thus, the shear correlations are isotropic. As a 
consequence of the relations ( |C.15| l, 

^±(i?)= d§'§'G^(-&,§')^,{^') (C.18) 
Jo 

with the survey field's 'radius' 6p ^ ^fA^Jn and the kernels (e.g. Crittenden et al. 20021 

G+(j?,»?') = \-6y){»-§'\ and (C.19) 



r^2n 1 r*2n r^oo 

G-(§,&')~ d^f>'G_(??Mo,)?'H^O = -r y e"^'''™"'^'^"cos[4^(0] d^p'e-'l'l'*' ^^I-^f^'-^'U dnum)Jo(M') 

Jo ^F ^ Jo Jo 

1 /4 12i?'^\ 

= ^M»?-t?') + (^-^)//('?-n 

where J denotes the Bessel function of the first kind, H the Heaviside step function, and 6d the Dirac delta function. 

Combining the relations ( |C.18| l between the convergence correlation and the shear correlation functions ^± with similar 
relations between and the correlations of the shear components y,-, one may derive (e.g. Kaiser|19"92 Schneider et al.|2002a| l 

UW) = {7i(O)yi(0 + + {y2(0)y2(0 + ^)} , (C.21a) 

= cos[4^{0)][{ym7i(O + &)} - (72(0)72(0 + t^))] + 2 sin[4^(0)] (n (0)72(0 + , (C.21b) 

(7dO,z)yi(0 + ^)} = ^-[UW)+cos[4<p(m-im), (C.21c) 

(72(0,2)72(0 + ^)} = ^^[Um - cos[4^(0)]^_(|,?|)), and (C.21d) 

<7i(0, 2)72(0 + sm[4<p(0)]UW). (C.21e) 



(C.20) 
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C.3. The shear correlation estimators 

We consider the following estimators |±(??) of the shear correlations ^±(§) at separations §: 

= with (C.22) 

Here, w*'^ denotes a statistical weight for the galaxy / (which we assume to be uncorrected with any other quantity of interest), and 
A{&, 0) denotes the bin window function ( [T4] i for a bin centred on ■& with width A{§) (which we assume to be small compared to 
the scale on which correlations change). Furthermore, y^'^' and yx ''' denote the tangential and cross component, respectively, of the 
shear y*'^ estimated from the shape of galaxy ; with respect to the direction = (f{6^^^ - 0*'*) defined by the positions (f'^ and d''-'^ 
of galaxies / and j. 

Since {j -i'^^ff^ + y^'-'^j'-f^},,, = ^±(|6'<-''' - 6'<''|) for / ^£ the estimator ( |C22l ) is an unbiased estimator of^J&) (e.g. 
let al.|2002b| l: 

\^''^^^/~\ 2;,.^w«wO)A()?,|eW-e»|) /~\ 2,jW»wO)A(i?,|0<^>-0«|) le'^^^^''' 

The sum Npi^) appearing in the denominator of the shear correlation estimators ( |C.22[ ) can be interpreted as the effective number 
of galaxy pairs in the bin centred on §. The expected effective number of pairs A^p(i?) in the bin reads 

Np(§) = {Np(§)) = V w«wW f ... f i^e^'^^* A(^, |0W - ^ ^A^p, (C.24) 
where 

AAi-^) ^ ^ f(fe f iVA(j?, \0' - e\) = 2n&A{d-) (C.25) 
Jf Jf 

denote the effective bin area and 

i.j=i 

denotes the effective total number of pairs. In the case of uniform weights w*^'^ - 1 V/, the effective number of pairs A^p = A^^, and 
the expectation ( |C.24[ ) reduces to: 

Np{i») = «^AfAa()?) = 2n§A(§)Ap n\. {C.21) 

C.4. The covariance of the shear correlation estimators 
We denote the covariance of the estimators ( |C.22| i by: 

c^±i»u^2) = {i±m,hm) - {U&i)){l±m) = i^^^^^^r^) - ^±(^1)^^21 (c.28) 

^ ' ^ ' \Npi&i)Np(^2)l 

To calculate c±±, we first neglect the randomness in the effective number of galaxy pairs Np{§y. 

A.p(.OA^p(^a) ^""-''^ 

The term 

{SMS.m) = I; w^'^^^^'w^''^'' (a(,?,, I^W - e(''|)A(^2, I^C) - (f% i^^'i^'-' ± (f?''>ff > ± f^'^rf ^)> (C.30) 



18 



Stefan Hilbert et al.: The log-normal approximation to cosmic shear covariances 

is then decomposed into three sums, exploiting the symmetry of the galaxy indices and the fact that only terms with even powers of 
intrinsic ellipticities contribute: 

(5.(^1)5,(^2)) = 2 1; (w<'>)' (a(,?i, - 0«|)A(^2, 100) _ 0(O|) (,(;V)^ao ^ ^a.)^ao) (^a;)^ao ^ ^a.o^ao)^ 

ij=i 

+ 4 ^ (wCf wWw«(a(^i , - 0® |)A(t?2, 10^'' - (e*!'^VT''''^ ± 'Vx"'') {ef'V^'^ ± V^'"'')) (C.3 1) 
ij,k=i 

The second sum comprises terms involving two or three distinct galaxies. The third sum contains terms involving two to four distinct 
galaxies. Since all galaxies are statistically indistuingishable apart from their weights w*'*, all averages {A(j?i, l^'-'^ - 0^'^\) . . .) in the 
same sum involving the same number of galaxies yield the same contribution to the covariance. However, terms with different 
numbers of galaxies involved yield different contributions. We neglected these differences, and thereby the noise due to the finite 
sampling of the shear field in position and redshift by the galaxies (which is justified for large A^g): 

{SMSM) - 2A^*^\A(,?„ 10(2) - 0")|)A(^2, 10(2) - 0(')|) (4''^^ef ') ± el!^\^^^'^) {ei^'^V^^^^ ± ei''\^^^'')) 

+4N,Ni^^\A{&, , 10(2) _ ^(i)|)A(^2, |e(^' - {erVr' ± 4''^?'^) (4''\?'" ± ei^'V^')) (C.32) 

+A^2(A(^, , 1^(2) _ ^(1)|)A(^,, 10(4) _ 0(3)|) (^0,2)^(2,1) ^ ^(1,2)^(2,1)) (^(3,4)^(4,3) ^ ^(3,4)^(43))^ 

The effective pair numbers 

Nf = ^ (^^'f Nf ^ = N-^' J] vvf^V® (C.33) 

i,j=l 

introduced here as abbreviations become identical to A^^ for uniform weights. 

The covariance obtained by combining Eqs. ( |C.29| l and ( |C.32| l is then split into three parts, with each part containing one of the 
three sums in the rh.s. of Eq. ( |C.32| l: 

c^±{'»u-»i) = cil(,?i,i?2) + cif (i?,,i?2) + c2'l('?i,t?2). (C.34) 
The ellipticity-noise part c"ll_ consists of all terms containing intrinsic ellipticities, but no shear: 

= ^^(^'^^^^^/ A(^i,|0(2) - 0(')|)A(^2, 1^^^) - {^^^'' ± 4''^'ef ')) (4'-\f'') ± e(''2)ef '))). (C.35) 

This ellipticity-noise part obviously vanishes for disjoint bins. The mixed part c"^]^ comprises all terms containing products of shear 
and intrinsic ellipticities: 

ctti^. ' '^2) - - «''l)^('^2, - {er\r' ± 4''2)^(2,l,) (41,3)^(3,1) ^ 41,3)^(3,1))^_ ^C.36) 

The cosmic variance part c^^ contains all terms involving only shear: 

c'lii^M = ^^(J;^^^^^ (A(>^i,|0(2> - 0(»|)A(^2, - e^^) (r<''^?'> ± r*''^?") (r?'\?'^' ± i^^i^"")) ^^.37) 

We proceed by evaluating the three covariance parts c^f) , c-^^ , and c'll for particular choices of the signs under various assumptions 
about the statistics of the convergence field. 

C.4.1. The ellipticity-noise contribution 



As mentioned above, c^fl vanishes for disjoint bins. The ellipticity-noise contribution to the covariance of ^±{d) evaluates to: 

c(f)(^,^) . ^^M^/(41.2),ai) ^ ^1,2)^(2,1)^41,2)^(2,1) ^ 41.2)^(2,m\ ^ 

A?p(l?)2 Ap \^ ^ ^ X X A T T XX)/ ^p(^)2 

Hence, the ellipticity-noise part also vanishes for non-matching signs. For uniform weights. 



f +(±l)(±l)f 



(C.38) 



c*,') (i?, )?) = c'fliff, §) = and c'lU^, §) = c''i{§, i?) = 0. (C.39) 



o", 



.4 
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C.4.2. The mixed contribution 

The mixed part ( |C.36[ ) comprises all terms containing products of shear and intrinsic ellipticities. After using the relations ( |C.21| i 
and some trigonometry, these reduce to: 

((,a,2)^(2,l) ^ ^(.,2)^(2,1)) (^(1,3)^(3,1) ^ 4''^<^.'))>^ = ^U\0^^^ - (C.40) 

+ 4'"V^^'")(4''\?'" - = Y^(l«*'' - O^'-'l) cos(4^<'-^> - 4^'^'^>), and (C.41) 

((^(1.2)^(2,1) _ ^,,,2)^(2,1)) (^(1,3)^(3,1) _ = ^^,(1^(3) - cos(4^('^^^) - (C.42) 

Hence, 



e o P 



3 r d^l' r d20(2) r d20«A(^,,|0P)_e(l)|)A(,?2,|0(3)_0(l)|)^^^(|0(3)_^(2)|) 

l,B Jf Jf Jf 2 

r d20(2) r d20(3)A(^^J0(2)|)^(^2,|0<^'|)^+(|0<^^-e<2>|) (C.43) 

Jf Jf 



AlNp{^i)Np{§2) 

r 

d(^l^+(|)?2Mo -)?lM^,|). 

Jo 



Here, denotes a unit vector with polar angle (p. Similar calculations yield: 



c^:i\&M = ^ ^ r d.^if_(|i?2«o - i?iH^J)cos[4^(»?2«^, - and (C.44) 



ci'^'(*^i,'?2)= ^ d^i^4|i?2«o-'?i«^,l)cos(4^i)- (C-45) 

A^n ^"g^F Jo 



'P 

For uniform weights, this reduces to: 

r.2 



c''tl\^u^i) = r d^if+(|i?2«o - (C.46) 

TTHgAp Jo 

c<f!'(i?i,»?2) = r d^i^_(|i?2«o - »?iH^J)cos[4^(i?2«^, - ^lU^M and (C.47) 

ci'^'(^i,»?2)= r d^i^+(|^2«o-t?i«^,l)cos(4^i). (C.48) 

nn^A^ Jo 

C.4.3. Tine cosmic variance contribution for general convergence fields 

To calculate the cosmic variance contribution to the covariance, we exploit the relations ( |C.15| l between the empirical correlation 
functions ( |C.12| l to transform the four-point correlations of the shear in the cosmic variance part ( |C.37[ ) into four-point correlations 
of the convergence: 

c'l^A^M--. ^ W^'^ fd'e*"' fd'e^'^ rd¥4)A(i^i,|0ra_0(i)|)A(,?2,|0'''-^<''l) 

AjAa(!?i)Aa(i?2) Jf Jf Jf Jf 

X i^y.^^'K 0(2) _ 0(I))^^(e(2), 0(2) _ 0(1)) + ^^(0(1), 0(2) _ 0(i))^,(0(2), 0(2) _ 0(1))] 

A^{§i)A^{i)-2) Jf Jf 

/± r l20(l)U (0(1), 0(2))^^(0(l) + 0(2), 0(2)) + ^^(0(1), 0(2))^^(0(l) + 0(2), 0(2))1 (C.49) 

\^F Jf 

d^ffC) fd¥2) fd^^* rd¥4)A(7?i,|0(2)|)A(,?2,|0(''l) 
3 Jf Jf Jf 



X I 



X 

'^F Jf 



AIA^(&i)Aa(§2) Jf 

fd20(2),(5^(^(2)^ ^(2),) r l2^(4),(5^(^(4)^ ^(4),) ^^(0(1))^(0(1) + 0(2)')^(0(3))^(0(3) + 0(4).)^^ _ ^^(^,)^^(^,). 
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From the properties of the kernels G± and G± follows that 

Wi) Jp Jf 



A/siifi) Jo Jo Jo Jo 

= -*- r d^<2) r*d#2)'^(2)' r V^'G^(^i„o,^(2)'„^p,)/(#2)'„^p„,...) (C.50) 
2:/!' Jo Jo Jo 

= — Pd#2)'^(2)'G^(^i,#2)') r'd.^P>^(^(2V^^^^, ) 

2;r Jo Jo 

= d^;^;G,(^i,^;)— — d0(2)'A(^;,|0(2)'|)/(0(2>,...) 

Jo '4a(e'\) Jf 



Thus, any contribution c^^ can be computed from a given C++ via a simple integral transform: 

Mr 



c. 



Jo Jo 



(C.52) 



C.4.4. The cosmic variance contribution for normal convergence fields 

For a normal convergence field, 

and thus (exploiting the symmetries of the integrands), 

c':liM2)=-^—^—— fd^e*" fd^^) fd¥3) rd20WA(^i,|0(2)|)A(^2,|0(4)|) 

A^a('9-iMa('?2) Jf Jf Jf Jf ^ ^ 

Jf Jf 

+ u\^^^ - eP\)U\^''^' + - - ^''D + + - - - - U^i)U^2) (C.53) 

AfAa(i?i)Aa()?2) Jf Jf Jf 

X r i20(2)'G^(0(2)^ ^(2V) r j2^(4V(5^(^(4)^ ^'^'o^.d^^'^ - e^^'>\)U\e^'> - e^"^\). 

Jf Jf 
Now, we make use of the following identity: 

r i^ff G-ie, ff)^^{W' -e'\)= \ iV G-{e, e" - 0')^Ki\o'\) 

Jf Jf 

Jf Af 

= cos[4^(6>) - 4^(6'" - 6»)] {d^ff e'-^^'-^-^^cos[A<p{O-A(p{0" -0)\^,{W\) (C-^^) 

Jf ^f 

[Aip{0) - Mff' - 0)] fdV y e'^(^"-^-«') sin[4.^(0 - 4.^(^' - 0)]UW\) 
Jf Af 



+ sm[ 

cos[4^c(6») - 4<piff' - 6»)]^_(|6»" - 0\). 



This can be generalized to: 

1 i^e' G±{G, ff)^^{W' - 0'\) = ^±{0, 0"), where (C.55) 
Jf 

0") = U\0" - 0\) and UO, 0") = ^-Q0" - 0\) cos[4ip(0) - 4ip(0" - 0)]. (C.56) 
Hence, 

AFAA(<7i)AA(f'2) Jf Jf Jf 



4 
^rAp 



* (C.57) 



Jo Jo 
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C.4.5. The cosmic variance contribution for zero-mean sliifted log-normal convergence fields 
For zero-mean shifted log-normal convergence fields, one can use Eq. ( |B.12| l to obtain 

+ [a\0^^^)U\0^'> - + 00) _ _ 0(2) _ 0(1)1)^,(1^(4) + 0(3) _ 0(2) _ 0(1)|) + . . .j 

+ K,'u\o'''m\(f'> - 0'''\)U\e''' + 0''' - 0'''\m0''' - 0''' - e^''\)U\0''' + 0^'' - 0'^' - 0'''\)UW''\). 

This reveals that the covariance for a zero-mean shifted log-normal convergence field equals the covariance for a normal convergence 
field (which consists of products of up to two correlation functions) plus 'corrections' involving products of three to six correlation 
functions. Inserting the coiTelation ( |C.58| l into the general expression ( |C.51| l for the cosmic variance contribution yields (again 
exploiting symmetries of the integrands) 



(C.58) 



C?](t?l,)?2) 



1 



A^Aa(i?i)Aa()?2) Jf 

X 



d20(i) rd20(2) rd20(3) [d^ef'^>^{§,,\e^^\)^{■&^^0(*^ 

' Jf Jf Jf 



I) 



i(3)r 



U\0^'^ - ^'\me''^ + 0(3) - 0(2) _ 0(I)|) + ^^(10(4) + 0(3) _ 0(1)|)^,(|0(3) _ 0(2) _ 0(1)|) + . . 

- , , fd'e^'' fd'0''' fd20(3)A(,?,,|0(2)|)A(^2,|0<'*)|){%„2^,(|0(2)|)^,(|0(4)|)^,(|0' 

ApAa(Vi)A^{§2) Jf Jf Jf >■ 

+ 2^„4[^2 ^ ^^(|0(2)|)j^^(|g(4)|^^^(|^(3)|)^^(|^(3) _ g(2)|) 

+ 2^o-4[^2 ^ ^^(|0(2)|)J[^2 ^ ^^(|0(4)|)J^^(|^(3) _ 0(2)|)^^(|0(3) _ 0(4)|) 

+ 4^„^[^2 + ^^(|0(2)|)][^2 ^ ^^(|0(4)|)j^^(|^(3)|)^^(|^(3) _ 0(2)|)^^(|0(3) _ 0(4)|) 

= 4KfU^l)U&2)^ Pd03 03^+(03) 
Jo 



(C.59) 



*ip, I) 



+ 2^01^0 +^+('^l)]^+('^2)^ r d03 03^+(03) r d^i^+(|03Ho-l?l«, 

Jo Jo 

+ 2^oV+('?l)[^0+^+('^2)]^ Pd03 03^+(03) r d^2^+(|03Mo-)?2«, 

Jo Jo 

+ 2^01^0 +^+('^l)][^0+f+('^2)]-r d03 03 d(^i^+(|03Ho-)?l«^,l) d^2 ^+(|03«O - 1?2«^J) 

^Ap Jo Jo Jo 

+ 4^o'^0+f+('^l)][^0+f+('?2)]4- r d03 03^+(03) r d^i^^l^jMo-t?!"^,!) f d^2 f+(|03«O " )?2«^J) 

^Ap Jo Jo Jo 

+ ^o'[^0+^+('^l)]k)+^+('^2)]z--r d0303^+(03) d^i^+(|03«O-»?l«^,l) d^2f+(|03«O-'?2«^,l) 

/ttAf Jo Jo Jo 

X ^+(|03«O - ^iU^i - ■&2U^,\). 

The expressions for c^ll and c'll can then be computed by combining Eqs. ( |C.51| l and ( |C.59| l. 

The expression ( |C.59|l f or a zero-mean shifted log-normal convergence field is considerably more complicated than the corre- 
sponding expression ( |C.53| l for a normal convergence field. If one only takes into account those terms in Eq. ( |C.59| ) that are also 
present in the case of a normal convergence field, and the simplest term oc k^^, one obtains 



rS-y) 



(.M2)^—r- d0303 dy>if±(|)?iH^,,03«ol) d(^2^±(|l?2«^„03«ol) + ^±(t?l)^±(»?2)^ d03 ^3 ^+(03)- 

^Ap Jo Jo Jo " /CqAf Jo 



(C.60) 



C.5. Shear correlation covariances for surveys comprised of several fields 

So far, we have discussed the case of a survey with a single rectangular field of view. If instead the survey consists of several 
independent fields of the same shape and size, and with the same galaxy density, the covariances for the whole survey can be 
obtained by (i) computing the covariances for a single field, and (ii) dividing the result by the number of fields. 
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The expressions derived here for the cosmic shear covariance depend on the survey geometry via the survey field area Ap. Its 
inverse A"' acts as an overall factor Furthermore, the linear dimension 9^ of the survey field enters as upper limit on integrals of 
correlation functions. Thus, the covariances for a survey do not only depend on the total survey area, but also on the shapes, sizes, 
and relative separatioiis of th e individual fields comprising the survey (in contrast to similar expressions derived by Schneider et al.| 



2002a 



Joachimi et al. 



2008 1. For example, the covariances for a survey with one 4 deg field usually differ from the covariances for 



a survey with 4 independent fields of 1 deg . 



Appendix D: Computing shear correlations with Fast Fourier Transform methods 

Our ray-tracing simulations provide us with the convergence and shear at positions forming a square mesh. Thus, Fast Fourier 
Transform (FFT) methods to compute their correlations are the obvious choice. Such FFT methods are so simple and efficient 
that they may be the best choice even for shear data sampled at irregular positions in fields with gaps and otherwise comphcated 
geometry (as in real shear surveys), where projection of the data onto a regular mesh is required as an intermediate step. 

As an illustrative example, we first revise how to use FFT methods to compute the two-point correlation of the projected galaxy 
density. Assume, the density contrast 6„{0) - {nJ(0) - (ng))/(«g), with rig denoting the projected galaxy density and (rig) denoting 
the mean projected galaxy density, is known on a subset of points 0, of a rectangular non-periodic mesh (0,), e.g. from a counts-in- 
cells. Some of the mesh points are allowed to have no data. To avoid aliasing, the mesh is first extended to a mesh twice as large as 
the original mesh in each dimension by adding mesh points around the original mesh. Let A^pix denote the number of mesh points of 
the padded mesh. 

An array w - (wi, . . . , wj^i^.J with a layout commensurate with the padded mesh geometry is used to store the window function 
w, which holds the information on which aiTay entries contain valid density data and what thek statistical weights (e.g. due to 
varying completeness) are. For example, 

f 1 if data for 0,- is vaUd, and , , 

[0 otherwise 

for uniform weights. Furthermore, the weighted density contrast data 

(w6g)i = w(0d6,(0d (D.2) 
on the mesh points 0i of the padded mesh are stored in an array wSg with the sam e geometry as the window array w 



An FFT algorithm (such as, e.g., provided by FFTW, Frigo & Johnson|2005 i is used to compute the discrete Fourier transform 
(DFT) w of the window array f] 



p 



The squared modulus of resulting array elements is then computed, 

= h.Aki) = w(kiTw(ki) = 2 e'*'-'''""''V(0^)w(0,), (D.4) 



m 

pix p,q 



and transformed back to configuration space to obtain the spatial auto-correlation S„^^ of the window function array: 
5 = S^AOi) = 7^ Z e'*' *''""''""^''w(0p)w(0,) = Z ^(0pX0^)^(0i^ 0P - (D.5) 



pix p,q,r p,q 



Here, 



fl if6» = 6»',and 

6(0,0') = . .' (D.6) 

[0 otherwise. 

The auto-correlation S„s,^wS, of the weighted density contrast is computed in the same way with FFTs. Finally, the correlations are 
binned into radial bins to obtain the text-book counts-in-cells Landy-&-Szalay estimator for the angular two-point coiTelation ^g^i 
of galaxies: 

. Zimm)s.s,..s,j ^ Zi.jm\0j-0i\M0iM0j)Ss(0i)Ss(0j) 

^""'^ ^" 2,. A(^,|0,|)5 " Zi,jm\0j-0i\M0dw(0j) 

where A{§, 9) denotes the bin function ([14]). 

With the same approach, one can calculate estimates for convergence and shear correlations. One creates arrays w, wk, wyi, and 
H'72 storing the values of the window function and the weighted convergence and shear components on the points 0/ of the padded 
mesh. One then uses FFTs to compute the auto-coiTelation S„ „ of the window function, the auto-coiTelation S^^.wk of the weighted 



The actual prefactor depends on the particular FFT implementation used. 
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convergence, and the auto- and cross-correlations S^yp.wy, of the weighted shear components. Finally, the correlations are binned 
into radial bins to estimate the isotropic correlations: 

, „ ^ i:imwi\)s..,..,i ^ i:,7 A(i?, \ej - g.i)wffl)w(g,-kfflMg,) 
^ i:iAi&,\Oi\)s„,.,i T.ijm\Oj-0i\>mw{ej) ' 

= 7- (D.9) 

~ i:ijm\oj-0i\)w(0dw(ej) ' 

Zi A(l?, I^il) {C0S[4^C(6',)]5 «,yi,wri,! - COS[4^C(0,)]5 wy2,wy2,i + 2 sin[4^(0,)]5 wyr,wy2,i} 

£-(■&) = : (D.IO) 

Y,iA(&,\ei\)s^AOi) 

_ 2,7 m \0j - OiDwiedwiOj) [yriOi, Oj - edy^iOj, 0j - 0i) - y^iOi, Oj - Ody^iOj, 0j - Od] 
~ j:ijAi&,\0j-Oi\)wi0i)wi0j) ■ 

FFT methods can also be used to compute galaxy-galaxy lensing statistics. For example, an estimator for the average tangential 
shear profile (yrK'S^) around a sample of lens galaxies as measured in a sample of source galaxies reads: 



'Zi A(l?, |6',|) COs[2^(0i)]5^,„,,„,,„,yi,,- -I- Sm[2(pi0i)]S^,„,,w,n,y2.i\ 
LiA{&,\0i\)Syt,^„^^v,ns,i 

where wi, n\, Ws, and denote the weights and densities of the lens and source galaxy samples. 
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